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THE SPIN-DOWN OF PSR J0821-4300 AND PSR J1210-5226: 
CONFIRMATION OF CENTRAL COMPACT OBJECTS AS ANTI-MAGNETARS 
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ABSTRACT 

Using XMM-Newton and Chandra, we measure period derivatives for the second and third known 
pulsars in the class of Central Compact Objects (CCOs) in supernova remnants, proving that these 
young neutron stars have exceptionally weak dipole magnetic field components. For the 112 ms 
PSR J0821-4300 in Puppis A, P = (9.28 ± 0.36) x ICT 18 . Its proper motion, ft = 61 ± 9 mas yr~\ 
was also measured using Chandra. This contributes a kinematic term to the period derivative via the 
Shklovskii effect, which is subtracted from P to derive dipole B s — 2.9 x 10 10 G, a value similar to 
that of first measured CCO PSR JI852+0040 in Kes 79, which has B s = 3.1 x 10 10 G. Antipodal 
surface hot spots with different temperatures and areas are deduced from the X-ray spectrum and 
pulse profiles. Paradoxically, such nonuniform surface temperature appears to require strong crustal 
magnetic fields, probably toroidal or quadrupolar components much stronger than the external dipole. 
A spectral feature, consisting of either an emission line at ~ 0.75 keV or absorption at s» 0.46 keV, 
is modulated in strength with the rotation. It may be due to a cyclotron process in a magnetic 
field on the surface that is slightly stronger than the dipole deduced from the spin-down. We also 
timed anew the 424 ms PSR J1210— 5226, resolving previous ambiguities about its spin-down rate. Its 
P = (2.22 ± 0.02) x 10~ 17 , corresponding to B s = 9.8 x I0 10 G. This is compatible with a cyclotron 
resonance interpretation of its prominent absorption line at 0.7 keV and harmonics. These results 
deepen the mystery of the origin and evolution of CCOs: why are their numerous descendants not 
evident? 

Subject headings: ISM: individual (Puppis A) ■ pulsars: individual (PSR J0821— 4300, 
PSR J1210-5226, PSR JI852+0040) — stars: neutron 



1. INTRODUCTION 

The class of faint X-ray sources in supernova remnants 
(SNRs) known as central compact objects (CCOs) are 
characterized by steady flux, predominantly surface ther- 
mal X-ray emission, lack of a surrounding pulsar wind 
nebula, and absence of detection at any other wave- 
length. Table [T] lists basic data on the well-studied 
CCOs, as well as proposed candidates whose qualifica- 
tions are not yet well established. Of the eight most 
secure CCOs, three are known to be neutron stars (NSs) 
with spin periods of 0.105, 0.424, and 0.112 s. Spin- 
down has been detected for two of these, the 0.105 s 
pulsar PSR J1852+0040 in Kes 79 and the 0.424 s pul- 
sar PSR J1210-5226 in the SNR PKS 1209-51/52. For 
PSR J1852+0040, t he implied surface dipole fie ld is only 
B s = 3.1 x 10 10 G (jHalpern fe GotthelJl2010al ). smaller 
than that of any other young known NS. In the case 
of PSR J1210— 5226, archival data allow two alternative 
timing solutions, with B s = 9.9 x 10 10 or 2.4 x 10 11 G 
HHalpern fc GotthelfllMl . 

It is natural to assume that CCOs that have not yet 
been seen to pulse are isolated, weakly magnetized NSs 
of the same class as the CCO pulsars. Where pulsar 
searches have been unsuccessful, it is possible that an 
even weaker magnetic field, a more uniform surface tem- 
perature, or an unfavorable viewing geometry, prevents 
detection of rotational modulation. The absence of pul- 
sations from the youngest known NS, the w 330 year 
old CCO in Cassiopeia A, has been used, in combina- 
tion with fitting of its X-ray spectrum, to argue that it 
is covered with a uniform temperature, non-magnetized 



atmospher e of carbon, the pro duct of nuclear burning of 
H and He (|Ho fc Heinkd 120091 ). Rapid cool ing of the NS 
inCas A, directly detected by Chandra (jHeinke fc Hd 
l2010f ). has been interpr eted as evidence for neutron su- 
perflu idity in the core (I Page et al.l 1201 It iShternin et al.l 

The "anti-magnetar" explanation of CCOs, which is 
motivated by their weak magnetic fields, absence of vari- 
ability, and location on the P — P diagram, remains 
incomplete in detail. Specifically, it does not yet ac- 
count for the hot spots that are seen on the surfaces 
of CCO pulsars. Since the spin-down power of a CCO 
pulsar is less than its X-ray luminosity, the latter must 
be thermal emission from residual cooling, which can 
only be nonuniform if there is anisotropic heat conduc- 
tion. In the absence of strong magnetic fields or mag- 
netospheric activity, it is difficult to reproduce the light 
curve and pulsed fraction of 64% from PSR J1852+0040 
in K es 79 (jHalpern fc GottheU I2010at iShabaltas fc Lail 
2012), or the two antipodal hot spots of different tem- 
peratures and areas on RX J0822— 4300 in Pu ppis A 
(|Gotthelf fc Halpernl [20091 IGotthelf et all [2010h . The 
latter 0.112 s pulsar, hereafter PSR J0821-4300, is a 
subject of this paper. Its spectrum is especially puz- 
zling in also displayi ng a phase-dependent em ission fea- 
ture at 0.7-0.8 keV (|Gotthelf fc Halpernl [2009h. which is 
repor ted to be variable in the long term (|De Luca et al.l 
I20i2h . 

Here we report the first spin-down measurement for 
PSR J0821— 4300, based on a dedicated program of 
phase-coherent timing jointly scheduled between XMM- 
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TABLE 1 

Central Compact Objects in Supernova Remnants 



ceo 


SNR 


Age 


d 


P 




B s 


L x ,bol 


References 






(kyr) 


(kpc) 


(s) 


(%) 


(10 iu G) 


(erg s i ) 




RX J0822.0-4300 


Puppis A 


4.5 


2.2 


0.112 


11 


2.9 


5.6 x 10 33 


1,2,3,4,5,6 


CXOU J085201.4-461753 


G266.1-1.2 


1 


1 




< 7 




2.5 x 10 32 


7,8,9,10,11 


IE 1207.4-5209 


PKS 1209-51/52 


7 


2.2 


0.424 


9 


9.8 


2.5 x 10 33 


6,12,13,14,15,16,17 


CXOU J160103.1-513353 


G330.2+1.0 


> 3 


5 




< 40 




1.5 x 10 33 


18,19 


1WGA J1713. 4-3949 


G347.3-0.5 


T.6 


1.3 




< 7 




~ 1 X 10 33 


11,20,21 


XMMU J172054.5-372652 


G350. 1-0.3 


0.9 


4.5 








3.9 x 10 33 


22,23 


CXOU J185238. 6+004020 


Kes 79 


7 


7 


0.105 


64 


3.1 


5.3 x 10 33 


24,25,26,27 


CXOU J232327.9+584842 


Cas A 


0.33 


3.4 




< 12 




4.7 x 10 33 


27,28,29,30,31,32,33 


2XMMi J115836. 1-623516 


G296.8-0.3 


10 


9.6 








1.1 x 10 33 


34 


XMMU J173203. 3-344518 


G353.6-0.7 


~ 27 


3.2 




< 9 




1.3 x 10 34 


35,36,37,38 


CXOU J181852.0-150213 


G15.9+0.2 


1-3 


(8.5) 








~ 1 x 10 33 


39 



Reference s. (1) IHui fc Becker 1120064 <2 \ JGotthelf &: Halpernl |2009|: (31 | Gotthelf et al.l 1201(3; (4) IDe Luca et alj|20ia (5) 

IBecker et "all 120121: (61 thi s paper: (71 Isiane et al-lpOOll: (81 IKargaltsev et al.l 120021: (91 IBamba et alJ 12003: (101 llvudin et al.M200a 
(11) Dc Luca 2008: (12) Zavlin ct al. 2000: (13) Mcrcghctti ct al. 2002; (14) Bienami ct al. 2003; (15) Dc Luca ct al. 2004; (16) 
Gotthelf & Halporn 2007; (17) Halporn & Gotthelf 2011; (18) Park ct al. 2006; (19) Park ct al. 2009; (20) Lazcndic ct al. 2003; (21) 
ICassam-Chenai et al.H2q04t ( 221 IGaensler et a.l.ll200St (2 31 ILoychinskv et al J 1 20111 ; (2 41 ISeward et al J 120031: (25 1 iGotthelf et alj|2005t 
(26l|Halpern et al.ll2007l: (271 IHalpern &: Gotthelfll2010at (2 81 IPavlov et a.l.ll2000t (2 91 IChakrabartv et al.l | 2001l: (301 IMereghetti et all 
2002; (31) Pavlov & Luna 2009; (32) Ho & Hcinkc 2009: (33) Hcinkc & Ho 2010; (34) Sanchez- Ayaso ct al. 2012; (35) Tian ct al. 2008; 
(36) Abramowski ct al. 2011; (37) Halpcrn & Gotthelf 2010b; (38) Halporn & Gotthelf 2010c; (39) Reynolds ct al. 2006. 

Note. — Above the line are eight well-established CCOs. Below the line arc three candidates. 

a Upper limits on pulsed fraction arc for a search down to P — 12 ms or smaller. 



Newton and Chandra. It was also necessary to incor- 
porate Chandra HRC observations of the position and 
proper motion of PSR J0821— 4300 in order to determine 
its small period derivative accurately. The astrometric 
analysis is described in Section 2. (The lat ter work was 
also c arried out, with consistent results, bv IBecker et al.l 
120121 ) The results of the timing are presented in Section 

3. The X-ray flux and spectra are discussed in Section 

4, with particular attention paid to the spectral line and 
the question of its possible variability. We also obtained 
new timing observations of PSR J1210— 5226 that resolve 
the prior ambiguity about its spin-down rate in favor of 
the smaller value; this definitive result is presented in 
Section 5. The nature of CCOs as anti-magnetars, and 
their possible evolutionary status, are discussed in Sec- 
tion 6. Conclusions and proposals for future work follow 
in Section 7. 

2. X-RAY POSITION AND PROPER MOTION 

Evidence that PSR J0821-4300 has high proper mo- 
tion from C handra HRC images ov er a 5 year baseline was 
report ed bv lHui fc Beckerl (|2006b| ) and (Winkler fc Petrel 
(2007), but with slightly disparate measurements of \x = 
107 ± 34 mas yr _1 and [i = 165 ± 25 mas yr^ 1 , respec- 
tively, from the same data. Here, we are concerned with 
timing this high-velocity pulsar over an extended period 
of time with millisecond accuracy. When trying to mea- 
sure a small P using X-ray data, position and proper 
motion can contribute significant errors via three effects. 
The first is an instrumental property of the Chandra 
CCDs when used in continuous-clocking (CC) mode (see 
Section 3); the position of the pulsar must be known a 
priori to < 0."5 in order to determine the time of arrival 
of each source photon. The second consideration is the 
accuracy of the barycentric correction. The third effect 
is the magnitude of the proper motion, which contributes 
a purely ki nematic period d erivative via the "train whis- 
tle" effect (|Shklovskiilll970D . The original measurements 



of proper motion are not accurate enough to measure this 
effect, which is crucial in the case of PSR J0821— 4300. 

Accordingly, we have reanalyzed the position and 
proper motion of PSR J0821— 4300 using the Chandra 
HRC-I data listed in Table [2j which now includes a more 
recent pair of observations in 2010 August that extends 
the baseline to 10.6 yr, enabling higher precision on both 
the contemporary position for timing, and the proper 
motion. We will describe here any differences between 
our method and previous work. For example, we did 
not use an HRC-S observation (ObsID 1851) because of 
known systematic differences between HRC-S and HRC- 
I. Ultimately, however, our results are consistent with the 
rec ent analysis of the s ame data (including ObsID 1851) 
bv IBecker etall (|2012t) . 

The data from all epochs were reprocessed and an- 
alyzed using the latest calibration files and software 
(CIAO 4.4/CALDB 4.4.8). This processing accounts 
for the H RC AMPJ5F electronic ringing distortions dis- 
cussed bv IHui fc Beckerl (|2006bD . The HRC detector is 
well suited to astrometry, with its processed pixel size of 
0."1318 that oversamples the on- axis point spread func- 
tion (PSF) by a factor of 5. For all observations, the 
pulsar was placed close to the optical axis where the PSF 
is essentially symmetric. In the following analysis we as- 
sume, as there is no evidence to the contrary, that the 
HRC focal plane is linear and the aspect reconstruction 
introduces no errors in roll angle. The two pointings on 
consecutive days in 2010 August are sufficiently differ- 
ent in their aspect reconstruction that we analyze them 
individually. 

The nominal uncertainty in aspect reconstruction for 
a typical Chandra observation is 0."6. It is often pos- 
sible to remove most of this systematic error by us- 
ing nearby X-ray point sources with precisely measured 
optical coordinates t o correct the absolute astrometry. 
IHui fc Bec ker (2006t3) used their X-ray detected "star A" 
as their sole fiducial point, and fitted a model PSF in- 
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TABLE 2 

Log of Chandra HRC-I Observations of PSR J0821-4300 



ObsID 


Date 


Start Epoch 


Exposure 


Roll angle 


Star A 


PSR J0821-4300 




(UT) 


(MJD) 


(ks) 


(°) 


(Counts) a 


(Counts) 3, 


749 


1999 Dec 21 


51533.95 


18.0 


338.7 


47 


3257 


4612 


2005 Apr 25 


53485.31 


40.2 


261.9 


123 


7260 


11819 


2010 Aug 10 


55418.72 


33.7 


163.4 


101 


5455 


12201 


2010 Aug 11 


55419.13 


38.7 


162.9 


117 


6296 



a Total counts collected in a l."5 radius aperture centered on the source. 



terpolated from the CIAO library appropriate for its po- 
sition a nd estimated pho t on en ergy to determine its po- 
sition. IWinkler fc Petrel (|2007l ) used this star and two 
additional stars, and followed the updated method out- 
lined in the CIAO thread for generating a Monte Carlo 
PSF using the CHaRT/MARX software for input into the 
CIAO/Sherpa spectral/image fitting software package. 

In our analysis, we also follow the CIAO thread to 
characterize the HRC-I PSF, but we adopt a simpler 
approach to measuring source locations, one that is not 
dependent on model fitting in the image domain. Our 
method is guided by the following observations. First, 
the on-axis, symmetric image of the pulsar contains 
enough counts that a simple centroid calculation is a 
sufficiently accurate measurem ent of its pos i tion o n 
the detector. Second, star A of iHui fc Becker! (|2006bh . 
which lies 2.'7 from PSR J0821-4300, is the only 
useful fiducial source for registering the X-ray image. 
The position and proper motion of star A are taken 
from the UCAC3 (jZacharias et al.l 120101 ). where it is 
listed as 3UC094-058669 with coordinates (J2000.0) 
R.A.=08 h 21 m 46. s 2924(16), decl.=-43°02'03."640(49), 
and proper motion fi a cos S = —14.3(2.0), fig = —3.6(5.5) 
mas yr . X-ray position m easurements of t h e two 
weaker, off-axis stars used by IWinkler fc Petrel ()20Q7f ) 
only add to the uncertainty (as quantified below) in the 
absolute astrometry. Third, the few X-ray photons from 
star A, and its broad off-axis PSF, do not warrant a 
sophisticated image fitting technique. Instead, we use a 
"corrected centroid" method, as described below. 

To determine the source location of star A in the X-ray 
images we start with the CHaRT/MARX simulation of the 
PSF, as described in the CIAO user webpagefl, for its re- 
spective locations on the focal plane. Figure [1] shows the 
distribution of the counts from each HRC image and cor- 
responding Monte Carlo PSF. It is apparent that star A 
is poorly sampled in the data, with total source counts 
in the range 47 < N < 123 (Table HJ. The maximum 
number of counts per pixel is typically only 2 — 4, mak- 
ing forward fitting poorly constrained statistically, while 
sharp features in the model can cause systematic offsets. 
Furthermore, the source is immersed in a substantial dif- 
fuse background from the Puppis A SNR which, although 
it only contributes a few photons over the source region, 
adds uncertainty to the position measurement because of 
the small count statistics, especially when fitting over a 
larger area. 

To sidestep these effects, the coordinates of star A are 
first found from its centroid. Photons were extracted 
from circular aperture of radius l."5, chosen to minimize 

1 http : / /exc . harvard . edu/chart 



the background counts. This aperture encloses essen- 
tially all of the signal in that fraction of the PSF with 
a finite probability of producing a single count during 
the observation. This measurement was made using the 
CIAO tool dmstat and is iterated to produce the final co- 
ordinates. However, while this results in a well-defined 
and statistically meaningful measurement, it does not ac- 
count for the shape of the complex off-axis PSF, whose 
orientation depends on the spacecraft roll angle (which 
differs for each observation; see Figure [1}, or for the 
monotonic off-axis deviations introduced by the flat focal 
plane. To quantify these systematic effects we also mea- 
sure the centroid of the simulated PSF of star A in sky 
coordinates in each image and compare it to the coordi- 
nates that were input to CHaRT/MARX for that PSF. The 
difference constitutes the small, but critical correction to 
the centroid of star A. 

We estimated the uncertainty in the derived coordi- 
nates of star A using a Monte Carlo method. We gen- 
erated 500 realizations of star A by sampling the PSF 
using a random number generator to match the observed 
counts, and accumulated the centroid measurements to 
build up a distribution in right ascension and declination. 
To account for the observed background, we included a 
random distribution of photons within the source aper- 
ture. The resulting (Gaussian) width of the distribution 
of centroids is typically a « 0."06 in each of the two co- 
ordinates (see Table |3|) . These reproduce the expected 
"standard error" for a centroid, s» a/y/N. We also simu- 
lated th e uncertainties fo r the t wo other fiducial sources 
used in IWinkler fc Petrel ()2007l ) and find that their in- 
clusion in the position determination would only increase 
the uncertainty on the final coordinates. 

The pulsar itself is a strong source (see Table whose 
coordinates are precisely measured using a centroid cal- 
culation with uncertainty an order of magnitude smaller 
than those measured for star A. Because of the symmetry 
of its on-axis PSF, no systematic correction is required 
for the pulsar. Its coordinates are then adjusted by the 
difference between the optical and X-ray coordinates of 
star A at each epoch to produce its final astrometric po- 
sition in Table [3] Fitting the position of the pulsar as a 
function of time then yields the proper motion. Figure [2] 
shows the \ 2 fit with constant velocity in each coordi- 
nate, and Table |4] lists the formal solution and quantities 
derived from it. The derived total proper motion of 61±9 
mas yr -1 is in agre ement with the value determined by 
IBecker et al.l mm . 71 ± 12 mas yr , and is smaller 
than previously published numbers for reasons discussed 
in that paper. 

The tangential velocity of PSR J0821-4300 then de- 
pends on a distance determination for Puppis A, which 
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Fig. 1. — Chandra HRC-I images around reference star A (3UC094-058669) used to define the coordinate system for position and proper 
motion of PSR J0821— 4300. Each panel shows the observed counts from Star A (left) and the simulated PSF (right) for its location on 
the focal plane in native HRC-I pixels. The plots cover 12" X 12" in celestial coordinates. The total counts in a l."5 radius aperture are 
given in Table [3] 



Spin-down of PSR J0821-4300 in Puppis A 



5 



TABLE 3 

Position Measurements for. PSR J0821— 4300 



Epoch 
(year) 


Star A (Optical) 
R.A. (h m s) Dccl. (° ' ") 


Star A (X-ray) 
R.A. (h m s) Decl. (° ' ") 


PSR J0821-4300 (corrected) 
R.A. (h m s) Decl. (° ' ") 


1999.975 
2005.316 
2010.610 
2010.611 


08 21 46.2928(16) 
08 21 46.2858(21) 
08 21 46.2789(31) 
08 21 46.2789(31) 


-43 02 03.637(49) 
-43 02 03.657(57) 
-43 02 03.676(76) 
-43 02 03.676(76) 


08 21 46.2882(65) 
08 21 46.3011(47) 
08 21 46.2621(62) 
08 21 46.2575(57) 


-43 02 03.397(83) 
-43 02 03.756(41) 
-43 02 04.113(47) 
-43 02 04.068(51) 


08 21 57.4024(67) -43 00 16.894(96) 
08 21 57.3685(52) -43 00 17.023(71) 
08 21 57.3488(69) -43 00 17.185(89) 
08 21 57.3476(65) -43 00 17.190(92) 



Note. — All coordinates arc equinox J2000. Optical coordinates for star A (3UC094-058669) arc corrected for the epoch of proper motion. 
The X-ray position of star A is determined using the method described in the text. The pulsar coordinates are corrected by the difference 
between the optical and X-ray coordinates of star A. Uncertainties on the last digits are in parentheses. 
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Fig. 2. — Four position measurements of PSR J0821— 4300 span- 
ning 10.6 yr, after correction using the optical/X-ray star A 
(3UC094-058669) as a reference. The two observations in 2010 
nearly coincide in time. The positions (diamonds) and their errors 
are fitted with a linear model (dashed lines) . The fitted parameters 
are listed in Table l4l 

has ranged from 1 to 2.5 kpc according to various meth- 
ods, and is a matter of unresolve d debate in the most 
recen t studies quoted here. One (jRevnoso et al.l H995l 
2003) is based on 21 cm H I velocities, and a perceived 
morphological association of features in H I surround- 
ing the pulsar and the SNR at v\ SI = +16 km s _1 , 
the latter corresponding to d = 2.2 ± 0.3 kpc. The 
secon d method uses spectra of ground-state hydroxyl 
lines (|Woermann et al.ll2000[ ). which show absorption at 
v isr < +7.6 km s" 1 , and emission above this velocity, 
from which d — 1.3+g'g kpc is derived. Both sets of 
authors employ assumptions that are not mutually ac- 
cepted, and are beyond the scope of this paper to evalu- 
ate. We will adopt a fiducial distance of 2.2 ± 0.3 kpc for 
purposes of further calculations, while noting the impli- 
cations of a possible smaller distance where relevant. 

The tangential velocity of PSR J0821-4300 at d = 
2.2 ± 0.3 kpc is vj_ = 636 ± 126 km s -1 . This velocity is 
at the high end of the distributi on of two-dimension al ve- 
locities of pulsars measured bv lHobbs et al. (2005), who 
find mean values of v± — 246 ±22 km s _1 for 121 ordi- 
nary (non-recycled) pulsars, and v± — 307 + 47 km s _1 
for 46 pulsars whose characteristic ages are < 3 Myr. 
The individual pulsar velocities are corrected for the so- 
lar motion with respect to the local standard of rest and 
a flat Galactic rotation curve in order to express them 
in the frame of the rotating Galactic disk. In the case 
of PSR J0821— 4300 this correction is dominated by the 
solar motion for the range of plausible distances, and is 
only —7 km s _1 , resulting in a corrected tangential ve- 
locity of v± yC — 629 km s _1 . If the distance is as small 
as 1 k pc, this reduces to an unexceptional 290 km s _1 . 

The Shklovskii (1970) effect, a purely kinematic con- 
tribution to the observed period derivative, will be sig- 
nificant. For a source moving at constant velocity the 
kinematic contribution is 



Pk = 



vjP 
dc 



(1) 



of 



From the proper motion measurement 
PSR J0821-4300, P k = (2.24 ± 0.72) x 10" 18 is 
calculated, where we have propagated the uncertainties 



TABLE 4 
Ephemeris of PSR J0821-4300 



Parameter 


Value 


Position and Proper Motion 


Epoch of position and /i (MJD) 


53964.0 


R.A. (J2000) 


08 h 21 m 57. s 3653(31) 


Decl. (J2000) 


-43°00'17."074(43) 


R.A. proper motion, /j, a cos8 


-54.1 ± 8.3 mas yr" 1 


Decl. proper motion, fig 


-28.1 ± 10.5 mas yr" 1 


Total proper motion, /i 


61.0 ± 8.8 mas yr" 1 


Position angle of proper motion 


242°5±9.°5 


Tangential velocity a , v± c 


629 ± 126 km s" 1 



Timing Solution 



Epoch of ephemeris (MJD TDB) b 
Span of ephemeris (MJD) 
Frequency, / 
Frequency derivative, / 
Period, P 

Period derivative, P 
Kinematic period derivative 11 , P^ 
Intrinsic period derivative 51 , P; nt 
Surface dipole magnetic field, B s 
Spin-down luminosity, E 
Characteristic age, r c 



55580.0000006 
55,182-56,027 
8.86529105448(32) Hz 
(-7.29 ±0.28) x 10" 16 Hz s" 1 
0.1127994550720(41) s 
(9.28 ± 0.36) x 10~ 18 
(2.24 ±0.72) x 10~ 18 
(7.04 ±0.80) x 10- 1S 
2.9 x 10 10 G 
1.9 x 10 32 erg s" 1 
254 Myr 



a Assuming d = 2.2 ± 0.3 kpc jRevnoso et al.N1995T ). 

k Epoch of fitted minima of the 1.5 — 4.5 keV pulse profile; phase zero in 
Figurc[5] 



on both fi and d. 



3. X-RAY TIMING 



Previous observations of PSR J0821— 4300 were 
only able to set upper l imits on its period d eriva- 
tive (IGotthelf & Halpernl [2001 IGotthelf et all |20T(1 
iDe Luca et all [2012] ) . Evidently P is so small that it 
can only be measured by phase-coherent timing. Ac- 
cordingly, we designed a sequence of observations coor- 
dinated between XMM-Newton and Chandra that would 
start and maintain phase connection over a 2 year span, 
2010 May - 2012 April. The scheduling strategy is 
the same as was used for P SR J1852+0040 in Kes 79 
(|Halpern fc Gotthelj l2010al ). By design, the resulting 
ephemeris was also connected backward to archival ob- 
servati ons that were obtaine d in 2009 December and 2010 
April (|De Luca et al.l 120121 ) , which extended the time 
span to 2.3 years. All of the timing observations used 
in this analysis are lis ted in Table El The disco very ob- 
servations from 2001 (jGotthelf fc Ha lpern 200J)|) are not 
included, as they are too far removed in time to be reli- 
ably connected. 

The Chandra observations used the Advanced Camera 
for Imaging and Spectroscopy (ACIS-S3) in continuous- 
clocking (CC) mode to provide time resolution of 
2.85 ms. All data were reprocessed from the level 1 event 
files with the coordinates corrected for the proper mo- 
tion of PSR J0821-4300 given in Tabled and analyzed 
using the latest calibration files and software (CIAO 
4.4/CALDB 4.4.8). Reprocessing with a source posi- 
tion that is accurate to < 1 pixel (< 0."5) ensures that 
the time assignment is precise to < 3 ms. All of the 
XMM-Newton observations used the pn detector of the 
European Photon Imaging Camera (EPIC-pn) in "small 
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TABLE 5 

Log of X-ray Timing Observations of PSR J0821— 4300 



Mission Instr/Mode 



ObsID Date Elapsed time/ Start Epoch Period a Z\ 

(UT) Livetime (ks) (MJD) (s) 



XMM 

XMM 

XMM 

Chandra 

XMM 

XMM 

XMM 

XMM 

XMM 

XMM 

Chandra 

XMM 

XMM 

Chandra 

XMM 

Chandra 

XMM 



EPIC-pn/SW 
EPIC-pn/SW 
EPIC-pn/SW 
ACIS-S3/CC 
EPIC-pn/SW 
EPIC-pn/SW 
EPIC-pn/SW 
EPIC-pn/SW 
EPIC-pn/SW 
EPIC-pn/SW 
ACIS-S3/CC 
EPIC-pn/SW 
EPIC-pn/SW 
ACIS-S3/CC 
EPIC-pn/SW 
ACIS-S3/CC 
EPIC-pn/SW 



0606280101 
0606280201 
0650220201 

12108 
0650220901 
0650221001 
0650221101 
0650221201 
0650221301 
0650221401 

12109 
0650221501 
0657600101 

12541 
0657600201 
12542,14395 
0657600301 



2009 
2010 
2010 
2010 
2010 
2010 
2010 
2010 
2010 
2010 
2011 
2011 
2011 
2011 
2011 
2012 
2012 



Dec 17,18 
Apr 05 
May 02 
Aug 16 
Oct 15 
Oct 15 
Oct 19 
Oct 25 
Nov 12 
Dec 20 
Feb 04 
Apr 12 
May 18 
Aug 11 
Nov 08 
Feb 18,19 
Apr 10 



85.1/54.9 
42.2/29.4 
28.0/19.6 
34.0/34.0 
23.5/16.4 
23.5/16.4 
26.5/18.6 
24.5/17.2 
23.5/16.5 
27.2/19.0 
33.0/33.0 
30.0/21.0 
36.5/25.6 
33.0/33.0 
37.2/26.1 
33.1/33.1 
35.3/24.7 



55182.820 
55291.377 
55318.782 
55424.625 
55484.109 
55484.987 
55488.332 
55494.228 
55512.524 
55550.159 
55596.837 
55663.857 
55699.925 
55784.655 
55873.289 
55975.446 
56027.022 



0.112799488(12) 
0.112799451(20) 
0.112799390(41) 
0.112799470(21) 
0.112799519(44) 
0.112799462(39) 
0.112799518(40) 
0.112799486(35) 
0.112799391(52) 
0.112799450(35) 
0.112799445(27) 
0.112799449(30) 
0.112799480(17) 
0.112799412(22) 
0.112799459(28) 
0.112799481(35) 
0.112799445(17) 



173.0 
99.1 
35.5 
92.3 
47.0 
56.7 
50.5 
62.1 
44.7 
63.3 
73.4 
56.6 
95.6 
80.4 
42.2 
65.8 
96.5 



Barycentric period derived from a Z\ test. The I Leahy et all dl983H uncertainty on the last digits is in parentheses. 



window" (SW) mode to achieve 5.7 ms time resolution, 
and an absolute uncertainty of rs 3 ms on the arrival 
time of any photon. Each data set was examined and 
cleaned of intervals of high particle background due to 
solar activity, as necessary. Two pairs of data sets that 
were acquired on consecutive days were merged to im- 
prove their statistics. The photon arrival times from all 
data were transformed to Barycentric Dynamical Time 
(TDB) using the coordinates of the pulsar corrected for 
proper motion. 

For PSR J0821-4300, the pulsed signal strength is a 
strong function of energy, not only because of its spec- 
trum relative to the background but because of the can- 
cel lation by emissi o n from the opposite pole, as described 
in iGotthelf et al.l (|2010T ) . For each observation listed 
in Table [5] we extracted source photons using an aper- 
ture centered on the source and optimized for the signal 
strength in the hard 1.5 — 4.5 keV energy band. For the 
XMM-Newton observations we used an aperture of ra- 
dius of 30". For the Chandra CC-mode observations we 
selected five columns (2."4). We also examined the soft, 
phase shifted 0.5 — 1.0 keV band. However, these data 
are noisier and their use did not in the end improve the 
timing results significantly. The phase cancellation ef- 
fect also prevents pulsations from being detected by the 
Chandra HRC, which has insufficient energy resolution. 

As in our previous t i ming studies of CCOs 
(jHafpern fc Gottnelfl l2010al 1201 ID . we employed two 
complementary approaches to fitting a n ephemeris . 
First, we used the Z\ (Rayleigh) test (|Struttl 118801 : 
IBuccheri et al.l 119831 ) in a coherent analysis of the entire 
set of 17 observations. Beginning with the closely spaced 
set spanning 2010 October 15-19, the Z\ test determined 
the pulse frequency with sufficient accuracy to connect in 
phase uniquely to the next observation. This procedure 
was iterated by adding each subsequent observation, and 
including a frequency derivative when it became evident. 
We also worked backward in time, incorporating all 17 
observations in the resulting unique ephemeris. The fit- 
ted frequency derivative is / = (— 6.94 ± 0.28) x 10~ 16 
Hz s -1 , where the lc uncertainty comes from the AZf = 



—2.3 contour around the peak power in (/, /) space. 

The second method also started with the Z\ test 
statistic, this time to find the period and pulse profile 
separately at each epoch. The 17 profiles were cross- 
correlated, shifted, and summed to create a master pulse 
profile template. The process was iterated to generate 
a more accurate template and a set of time-of- arrival 
(TOA) measurements and their uncertainties for each 
epoch. These TOAs were fitted with a quadratic model 
in frequency and frequency derivative using a \ 2 fitting 
routine to minimize their phase residuals. We searched 
for a coherent phase-connected solution over a grid of / 
and / covering the range / = 8.8652906 ± 0.0000016 Hz 
(at epoch MJD 55,580) and -3.1 x 10~ 14 < / < 1.9 x 
10~ 14 Hz s _1 , with an oversampling factor of 10 for ac- 
curacy. This range corresponds to the 3cr limits of an in- 
coherent fit to all of the measured frequencies, including 
the 2001 discovery observations. The resulting frequency 
derivative from TOA fitting, / = (-7.29 ± 0.28) x 10~ 16 
Hz s _1 , is consistent with the value found above from 
the coherent Z\ search. We adopt the TOA result for 
the final timing solution listed in Table 2J 

The observed P = (9.28±0.36) x 10~ 18 can now be split 
into the sum of its intrinsic and kinematic contributions, 
P = Pi n t + Pfe. Since we determined in Section 2 that 
P fc = (2.24±0.72) x 10~ 18 , the intrinsic period derivative 
is P int = (7.04 ±0.80) x 10" 18 . Parenthetically, we note 
that the small observed period derivative is independent 
evidence that the proper motion of th e pulsar is not as 
large as the value originally quoted bv lWinkler fc Petri 
(2007), ii = 165 ± 25 mas yr" 1 . If so, and if d = 2.2 kpc, 
Pfe would be (1.64 ± 0.55) x 10~ 17 , requiring P int to be 
negative, i.e., the pulsar would be spinning up. 

In the vacuum dipole spin-down formalism, the values 
of P and Pnt imply a surface magnetic field strength 
B s = 3.2 x lO^PP) 1 / 2 G = 2.9 x 10 10 G, a spin- 
down luminosity E = -Jfifi = 4tt 2 PP/P 3 = 1.9 x 

10 32 erg s _1 , and characteristic age t c = P/2P = 
254 Myr. PSR J0821-4300 is nearly identical in its spin 
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Fig. 3. — P — P diagram of isolated pulsars (dots), binary radio 
pulsars (circled dots), and other types of isolated X-ray pulsars 
(colored symbols). The CCO pulsars (red stars) in Kes 79 and 
Puppis A have virtually the same spin parameters. The upper 
limit on Calvera's P is from Halpcrn (2011). The radio pulsar 
death line B/P 2 = 1.7 X 10 11 G s~ 2 oflBhattacharva et ail 1T992T) 
is indicated. The van den Heuvell (fl9S7(T spin-up limit for recycled 
pulsars corresponds to P(ms) = 1.9(S/10 9 G) 6 / 7 . The exponent 
in this equ ation corrects a typ o graphi cal error in the caption to 
Figure 7 of Halpern & Gotthelf (2010a), although the correspond- 
ing line in the Figure was correct. 



properties to PSR J1852+0040 in Kes 79, as shown in 
Figure [3] The uncertainties in distance and proper mo- 
tion have only a small effect on the derived magnetic 
field. For a smaller distance of 1 kpc, is reduced 
to 1.02 x lCT 17 , and B s = 3.1 x 10 10 G. An absolute 
upper limit regardless of distance and proper motion is 
B s < 3.3 x 10 10 G. 

The phase residuals from the best-fit solution are 
shown in Figure |4] The weighted rms of the phase residu- 
als is 5.1 ms, or 0.045 pulse cycles, which is comparable to 
the individual measurement errors (average a = 3.6 ms). 
It is not clear if there is any real timing noise and/or 
systematic errors in the TOAs. 

The light curves of PSR J0821-4300 in the soft and 
hard bands are shown in Figure [5j These were derived by 
folding all the timing data on the best fitting ephemeris 
given in Table |H As revealed by a cross-correlation, the 
soft and hard pulses are out o f phase by 0.45 ±0.02 cycles, 
consistent with that found bv lDe Luca et all (|2012[ ) using 
the 2009 December and 2010 April XMM-Newton data. 

The energy dependence of the pulsar modulation pro- 
vided an important diagnostic for modeling the view- 
ing geometry and sur face emission of PSR J0821— 4300 
(|Gotthelf et al.l I2010T) . By analyzing the lightcurve in 
narrower energy bands than in Figure [5J we can resolve 
the signal modulation and phase for PSR J0821— 4300 
over the 0.3—5 keV range. The data were grouped into 23 
energy bands that are at least 1 00 eV in width and have 
a signal-to-noise Ng/^Ng + N b > 100, where N s ,N b are 
the source and background counts, respectively. We used 
Z\ to provide a model of the unbinned lightcurve. The 
first Fourier component is a reasonable estimate as the 
lightcurve is sinusoidal in each energy band to within 
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Fig. 4.— For the timing observations of PSR J0821-4300 listed 
in Table \E\ pulse-phase residuals from the linear term (dash-dot 
line) of the phase ephemeris presented in Table [3] The quadratic 
term (solid line) contributes ss ±0.5 cycles to the ephemeris over 
the 2.3 years of timing. 

counting statistics. The error bar for the phase is cal- 
culated by cross-correlation, with the profile of Figure [5] 
serving as a template. 

The result is presented in Figure [6j As the energy 
dependent modulation decreases, the phase becomes un- 
defined in two energy bands; these two phase points are 
not plotted. The modulation is qualitati vely similar to 
that p redicted by the antipodal model of iGotthelf et al.l 
(2010) (cf. their Figure 6), providing confirmation of 
the basic model. The prediction was based on fitting 
the modulation in only three energy bands, using much 
less data, and differs somewhat from the new, resolved 
data having an order-of-magnitude more counts. In par- 
ticular, the energy of the minimum modulation is lower 
(1.12 vs. 1.28 keV), and the form of the modulation is 
more complex than predicted. Furthermore, the phase is 
seen to drift at lower energies and the transition is not as 
sharp compared to the antipodal case, in which the phase 
was statistically either 0.0 or 0.5. The observed charac- 
teristics likely imply that the hotspots are offset from a 
strictly antipodal geometry. The high quality data pre- 
sented herein should allow a far more detailed modeling 
of the surface emission of PSR J0821-4300. 

4. SPECTRAL ANALYSIS 

Previously, we modelled the original (year 2001) 
XMM-Newton spectra of PSR J0821-4300 as surface 
blackbody emission from two antipodal spots of differ- 
ent temperatures and areas. Crucially, this model is also 
able to account for the observ ed energy-dependent pulse 
modulation and pha se shift (|Gotthelf fc Halpernl 120091 : 
IGotthelf et aD l2010h . Additionally, we found evidence 
of a spectral line featur e around 0.8 k e V, w hich more 
recent data obtained by iDe Luca et al.l (|2012f) suggests 
has a time- variable centroid energy. With the increased 
quantity of data now in hand on PSR J0821— 4300, we 
can re-examine this spectral feature, first by combining 
all 13 XMM-Newton observations presented in Table [5] 
plus the two observations obtained in 2001, and then by 
testing for any variability among the observations. 

4.1. Summed XMM-Newton Spectrum 
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Fig. 5.— Summed pulse profiles of PSR J0821-4300 in the 1.5 - 
4.5 keV band (top) and the 0.5 — 1.0 keV band (bottom) using 
all of the observations listed in Table \E\ folded according to the 
ephemeris of Table [4] These hard and soft pulse profiles are out of 
phase by <j> = 0.45 zt 0.02 cycles. The intervals between the vertical 
lines (A0 = 0.3 cycles) correspond to the two phase regions used 
in the phase resolved analysis of Section 4.3 and Table [3] 
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Fig. 6. — Modulation and pulse phase as a function of energy 
for PSR J0821-4300 using all of the observations listed in Tabled 
folded according to the ephemeris of Table [31 The modulation 
(top) and pulse phase (bottom) reproduce the form predicted for 
the antipodal model of Gotthclf ct al. (2010). Two undefined phase 
points are omitted. 

The spectral analysis presented here uses exclusively 
the data collected with the EPIC pn. Although back- 
ground is much reduced in the EPIC MOS, this instru- 
ment is less sensitive at soft energies, collecting 4.7 times 
fewer photons than the EPIC pn in the 0.3 — 1 keV band. 
Furthermore, much of the EPIC MOS data were lost as 
the high surface brightness of the Puppis A SNR fre- 
quently triggered an automatic shutoff of that detector. 
We also neglect the Chandra ACIS spectra here because 
of their even poorer low-energy sensitivity and the in- 
creased background and other uncertainties involved in 
the analyzing the Chandra CC-mode data. Lastly, one 
ACIS image taken in timed exposure mode (ObsID 750) 
suffers from pileup, and was not used. 

For the EPIC pn data, the main technical issue is that 
the SNR background exceeds the point source counts be- 
low 1 keV. Surprisingly, the background intensity and 
spectral shape are both strong functions of spacecraft 



roll angle. This effect is made evident because the XMM- 
Newton data sets were acquired in two narrow ranges of 
roll angle roughly 180° apart associated with their re- 
spective visibility windows, with the time divided nearly 
equally between the two. We checked all of our spectral 
results carefully for features that might be dependent on 
roll angle due to systematic errors in background sub- 
traction. No such systematic effect were found, which 
indicates that the background subtraction is reliable in 
general. 

We extracted spectra for each observation from the 
EPIC pn detector using an aperture of radius Of 3 and a 
concentric background annulus of 0.'5 < r < 0.'6, select- 
ing only events with PATTERN < 4 and FLAG= 0. The 
data were filtered to exclude time intervals of high back- 
ground identified by count rate > 0.1 s _1 in the 10—12 
keV energy band. An inspection of the pattern distri- 
bution of single and double events shows no evidence 
of pile-up and suggests that a lower energy bound in the 
range 0.3 — 0.4 keV is acceptable. Response matrices and 
effective area files were generated for each observation 
using the SAS software suite. We combined data from 
all observations using the FTOOL addascaspec to pro- 
duce a single source spectrum and associated files. The 
combined spectrum was grouped to include at least 1000 
counts per channel and was fitted using XSPEC vl2.21 to 
a two blackbody model with interstellar absorption over 
the energy range 0.3 — 5 keV (see Figure [7] and Table [6]). 

The resulting high signal-to-noise ratio of the fit- 
ted spectrum reveals significant features that are unac- 
counted for by the two blackbody model. The best fit 
model, with N H = (3.8 ± 0.01) x 10 21 cm" 2 , kT w = 
0.26 ± 0.01 keV, and kT h = 0.46 ± 0.01 keV, has reduced 
xl = 1-50 for 359 degrees of freedom (DoF), which is for- 
mally unacceptable. The deviations are evident in struc- 
ture in the residuals in Figure [7K- Adding a Gaussian 
emission line to the model as suggested by our previous 
work improves the fit to xt = 1-08 for 358 DoF (Fig- 
ure^). The centroid energy of the line is 0.75±0.01 keV, 
and its equivalent width is 53 ± 10 eV. 

Considering the shape of the residuals from the two 
blackbody fit in Figure [7^, an alternative hypothesis 
is that an absorption feature, or features, is respon- 
sible. Accordingly we applied a Gaussian absorption 
line, two Gaussian absorptio n lines, and final l y, the 
cycl otron absorpt i on mo del of [Makishi ma et al.l (1990) 
and iMihara et all (|1990l ). which is available in XSPEC as 
cyclabs. As shown in Tableland Figure all of these 
models gave acceptable fits. Two absorption lines, when 
fitted independently, are separated by nearly a factor of 
2 in energy, which is suggestive of the fundamental and 
first harmonic in a cyclotron model. Therefore we fixed 
the ratio of their centroid energies at 2, with the result 
that their centroids are at 0.46 keV and 0.92 keV, brack- 
eting the energy previously ascribed to an emission line. 
The cyclabs model, which also fixes this ratio, yields 
the same centroid energies. We conclude that, from the 
phase-averaged spectra alone, it is not possible to dis- 
tinguish an emission feature from one or two absorption 
lines, which leaves the physical interpretation uncertain. 

The fundamental energy of the electron cyclotron res- 
onance falls at 



E = 1.16 {B/l0 iL G)/(l + z) keV, 



(2) 
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TABLE 6 

Models for the Summed XMM-Newton Spectrum of PSR J0821-4300 



Model 




Two Blackbody 


+ Emis. line 


+ Abs. line 


+ Two Abs. lines 


+ Cyclabs 


7V H (10 21 cnT 




3.8 ±0.1 


4.3 ±0.3 


3.2 ±0.2 


2.9 ±0.4 




2.8 ±1.01 


kT w (keV) 


0.26 ± 0.01 


0.25 ± 0.01 


0.29 ± 0.01 


0.29 ±0.01 




0.28 ± 0.03 


kT h (keV) 




0.46 ± 0.01 


0.45 ± 0.01 


0.49 ± 0.02 


0.49 ± 0.03 




0.47 ±0.02 


L™(bol) (10 33 


erg s ) a 


3.3 ±0.1 


3.6 ±0.2 


3.1 ±0.2 


3.0 ±0.2 




3.2 ±0.8 


L h (bol) (10 33 


erg s _1 ) a 


2.0 ±0.2 


2.0 ±0.2 


1.4 ±0.3 


1.4 ±0.3 




1.7 ±0.4 


A m (km 2 ) 




72 ± 11 


89 ± 18 


44 ±6 


40 ±8 




53 ± 14 


A h (km 2 ) 




4.4 ±0.9 


4.5 ±0.9 


2.5 ±0.8 


2.4 ± 1.0 




3.4 ±0.8 


E (keV) 






0.75 ±0.01 


0.46 ± 0.05 


0.46 ±0.01 




0.46 ± 0.01 2E 


Width b (oV) 






75 ±20 


85 ±50 


106 ± 20 


34 - 62 


53 - 97 35 - 290 


EW (eV) 






53 ± 10 










T c 








0.1 - 0.9 


0.6 - 1.3 


< 0.035 


0.9-1.7 <0.14 


x 2 (DoF) 




1.50(359) 


1.08(356) 


1.16(356) 


1.08(354) 


1.12(354) 



Note. — The lcr uncertainties for three interesting parameters (Ax — 3.53) are given. 
a Blackbody bolomctric luminosity for a distance of 2.2 kpc. 

k Gaussian a for the emission or absorption lines, natural width W for the cyclotron absorption model (Makishima ct al. 1990; 
IMihara eraT||1990|) . 
c Optical depth at line center. 
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Fig. 7. — (a) EPIC pn spectrum of the 16 summed XMM- 
Newton observations of PSR J0821— 4300 fitted to a double black- 
body model. The residuals from the fit are shown in the bottom 
panel, (b) The same spectrum fitted with a double blackbody 
model plus Gaussian emission line. The parameters of these fits 
are given in Table [6] 



ENERGY (KEV) 

Fig. 8. — Same as Figure [7] but showing the residuals from 
fits to a double blackbody model plus, from top to bottom: no 
line, a Gaussian emission line, a Gaussian absorption line, two 
Gaussian absorption lines, and a cyclotron absorption line model. 
The parameters of the fits are given in Table [6] 

where z is the gravitational redshift. Assuming a typical 
value of z = 0.3, and if B « B s = 2.9 x 10 10 G, the 
equatorial held from the vacuum dipole spin-down result, 
then Eq sa 0.26 keV is expected. When compared with 
Eq = 0.46 keV from the absorption-line ht, or 0.75 keV 
from the emission-line fit, this hints that the local surface 
field where the line is formed is larger than the equatorial 
dipole field, but is perhaps more compatible with the field 
at the pole, B p = 2 B s . Other factors possibly affecting a 
comparison between the dipole spin-down _B-field and the 
surface B-ficld from the cyclotron energy will be explored 
in Section 6.2. 

The addition of emission or absorption lines to the 
model near the low-energy end of the XMM-Newton 
spectrum affects the fitted value of the column density 
TVh, with values ranging from (2.8 — 4.3) x 10 21 cm~ 2 for 
the different models in Table [51 Therefore, it is possi- 
ble that comparison with independent measurements of 
Nk may suggest a preference for one or another of these 
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Fig. 9. — Blackbody temperatures and flux in the 0.5 — 5 
keV band for the 16 individual XMM-Newton observations of 
PSR J0821-4300, fitted to the two blackbody model. Data are 
taken from Table[J] Errors bars are la. The weighted mean values 
are indicated by the dashed lines. 

models. For example, the 21 cm H I emission in the fore- 
ground of P uppis A amoun t s to T Vht = 2.5 x 10 21 cm" 2 
according to lRevnoso et al.l ()2003[ ). To obtain this value, 
they integrated the H I line emission in the radial velocity 
range —10 to +16 km s , the latter velocity correspond- 
ing to their assumed 2.2 kpc distance of Puppis A. Com- 
parison with the values of Ah in Table [5] tends to favor 
the X-ray spectral models that include absorption lines. 
While this agreement is encouraging, it is subject to the 
caveat that the 21 cm column density is formally a lower 
limit, assuming as it does that the line is optically thin. 
Better support for a low column density comes from the 
XMM-Newton RGS spectra of seve r al reg ions of the Pup- 
pis A SNR fitted bv lKatsuda et al.l (|2012f ). These require 
X-ray A H in the range (2.58 - 2.85) x 10 21 cm -2 , which 
also agrees closely with the Ah from our absorption-line 
models for the pulsar spectrum. 

Finally, we note that the total luminosity and black- 
body are as reported in Tabl e EH di ffer from those pre- 
sented in iGotthelf &: Halpernl (|2009f) , the areas by a fac- 
tor of two or more. This is a consequence of the best fit 
values for each model fit. The derived blackbody areas 
depend on blackbody temperature as T~ 4 , which itself is 
strongly correlated with the fitted column density. Fu- 
ture work, simultaneously modeling of the phase depen- 
dent spectra, should better constrain the column density 
and temperatures, and consequently provide a more ac- 
curate measurement of the blackbody areas. 

4.2. Search for Variability 

To test for long-term variability of PSR J0821-4300, 
we also fitted the individual XMM-Newton spectra, 
searching for temperature variations on the surface, for 
example. The spectra for each observation are well-fitted 
by the two blackbody model, with parameters listed in 
Table [7j their statistics do not warrant an additional line 
component. For this comparison, we held the absorbing 
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Fig. 10. — Combined XMM-Newton EPIC pn spectrum of 
PSR J0821— 4300 in two pulse-phase intervals defined in Figure [5] 
(A</>i, blue, A02i red), fitted to a double blackbody model plus 
Caussian emission line. The blackbody temperatures and line cen- 
troid energy are linked between the two spectra. The fitted param- 
eters are given in the last column of Table [8] The residuals to the 
fit are shown in the bottom panel. 

column fixed at the value in Table [5] derived from the fit 
of the composite spectrum to the two blackbody model. 
Figure [9] displays the results. No significant variation 
is evident in the flux or in either of the two blackbody 
temperatures. 

We next test for pulse-phase dependence of the spec- 
tral line feature(s) by generating spectra in two phase 
intervals, each of A(j> = 0.3 cycles in width centered on 
the peaks of the soft and hard light curves, respectively 
(see F igure E) This is motivate d by the original find- 
ing in IGotthelf fc Halpernl fj2009f ) that an emission line 
is more strongly associated with the warm region, i.e., 
the soft phase of the pulse. For simplicity, given the re- 
duced counts in the phased spectra, we used only the 
Gaussian emission-line model as a representative for all 
of the possible line models. The two phase-resolved spec- 
tra combined from all epochs are shown in Figure 1 1 01 
where they are fitted to the two blackbody plus Gaus- 
sian emission-line model. In the simultaneous fit, the two 
temperatures and the line centroid energy are linked. As 
shown in Tabled the equivalent width of the Gaussian 
line is indeed about factor of 2 larger in the soft-phase 
s pectrum than in the ha rd phase. 

iDe Luca et al.1 ((2012) presented evidence that the 
emission-line centroid had decreased from 0.80 keV to 
0.73 keV between 2001 and 2009. Here we extend the 
examination for variability of the spectral feature, in- 
cluding its phase dependence, by combining the XMM- 
Newton spectra into four groups of adjacent observations 
as indicated Table [7] We fitted these four sets of spectra 
to the two blackbody plus Gaussian emission-line model, 
with their temperatures and line width (sigma) linked 
between the two phase intervals. The four fits are dis- 
played in Figure [TTJ and the fitted parameters for each 
set and their sum are presented in Table [5J It is clear 
that the temperatures and their contributed fluxes are 
steady in time, consistent with their phase-averaged be- 
havior and with no time dependence in their phase ra- 
tios. On the other hand, the equivalent width of the 
fitted Gaussian emission line shows evidence for having 
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TABLE 7 

Individual XMM-Newton Spectra of PSR J0821- 



4300 



Group 


ObsID 


Livetime 
(ks) 


(keV) 


(keV) 


Flux (10~ 12 
erg s — 1 cm -2 ) 


X?(DoF) 


1 


0113020101S 
0113020301S 


15.2 
16.1 


0.27(0.25,0.29) 
0.27(0.25,0.28) 


0.47(0.44,0.52) 
0.46(0.43,0.50) 


4.23(4.14,4.30) 
4.21(4.12,4.27) 


1.06(154) 
1.05(184) 


2 


0606280101S 
0606280101U 
0606280201S 
0650220201S 


25.7 
24.5 
24.6 
8.1 


0.27(0.26,0.28) 
0.25(0.24,0.26) 
0.27(0.26,0.28) 
0.22(0.19,0.26) 


0.47(0.45,0.51) 
0.45(0.43,0.47) 
0.47(0.45,0.50) 
0.40(0.38,0.45) 


4.15(4.09,4.20) 
4.11(4.05,4.16) 
4.14(4.09,4.20) 
4.11(4.01,4.20) 


1.19(282) 
0.91(268) 
1.02(261) 
0.88(94) 



0650220901S 
0650221001S 
0650221101S 
0650221201S 
0650221301S 
0650221401S 



16.4 
15.5 
18.6 
17.1 
16.4 
15.5 



0.24(0.22,0.25) 
0.26(0.25,0.28) 
0.27(0.26,0.28) 
0.28(0.27,0.29) 
0.24(0.23,0.26) 



0.27(0. 
0.26(0. 
0.26(0. 
0.26(0. 



25,0.28) 
25,0.27) 
25,0.27) 
25,0.28) 



0.42(0 
0.48(0 
0.49(0 
0.52(0 
0.43(0 



0.45(0. 
0.47(0. 
0.45(0. 
0.46(0. 



42,0.< 
40, 0.44) 
45,0.51) 
46, 0.54) 
48,0.57) 
41,0.46) 



43,0.48) 
44,0.50) 

43, 0.47) 

44, 0.49) 



4.13(4.06,4.20) 
4.19(4.12,4.25) 
4.17(4.08,4.23) 
4.06(3.98,4.12) 
4.13(4.05,4.20) 



4.26(4. 
4.13(4. 
4.24(4. 
4.22(4. 



20, 4.32) 
07,4.19) 
19, 4.30) 
16,4.27) 



0.99(183) 
0.97(222) 
1.08(201) 
1.21(182) 
1.00(174) 



0.96(225) 
0.99(242) 
1.07(285) 
0.94(262) 



0650221501S 
0657600101S 
0657600201S 
0657600301S 



20.3 
22.4 
26.1 
24.7 



Note. — Results for a fit to an absorbed, two-blackbody model with the column density held fixed at 
Nh — 3.75 X 10 21 era -2 , the value obtained for the summed spectrum in Table[6] The absorbed flux in the 
0.5 — 5 keV range is tabulated. The range of uncertainty of each value is the la confidence interval for two 
interesting parameters. 




Fig. 11. — XMM-Newton EPIC pn spectra of PSR J0821-4300, grouped into the four sets as listed in Table0 and fitted to a double 
blackbody model with Gaussian emission line. As in Figure flOl red and blue represent the two pulse- phase intervals defined in Figure [5] 
The residuals from the fits are shown in the bottom panels. The fitted parameters are given in Table [8] 
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increased between 2001 and 2009, while the line centroid 
energy decreased from 0.79±£gl keV to 0-71 + ° ° 2 keV. 
While these results are consistent with the lDe Luca et al\ 
<|2012f ) analysis of the same data, there is little if any ad- 
ditional variability between 2009 and 2012. 

We also repeated this test for variability with the Gaus- 
sian absorption-line model for the spectral feature (re- 
sults not tabulated here) . For either the emission-line or 
the absorption-line model, the measured line centroids 
for all epochs but 2001 cluster well within their lcr un- 
certainty (for two interesting parameters, A^ 2 = 2.3). 
The 2001 set deviates from the mean defined by the other 
three sets by « 2a. In terms of percentage deviation from 
the mean, the line centroid measured in 2001 was 14% 
and 8% higher in the emission line and the absorption 
line model, respectively. We conclude that, although the 
deviation seems large for the emission line model, it is 
not inconsistent with the expected variance. Further ob- 
servations would be necessary to establish more definite 
evidence of long-term variability. 

5. A DEFINITIVE SPIN-DOWN MEASUREMENT 
FOR PSR J1210-5226 

Archival timing observations of PSR J1210— 5226 
spanning the years 2000-2008 were too sparse to se- 
curely determine it s spin -down rate, as described in 
lHalpern fc GottheH ()2011[ ). Searching all possible pa- 
rameter space for a phase-coherent, quadratic ephemeris, 
we found two equally acceptable solutions, with / = 
-1.243(22) x 10~ 16 Hz s" 1 and / = -7.084(22) x 10~ 16 
Hz s _1 , corresponding to B s = 9.9 x 10 10 G (solution 1) 
and B s = 2.4 x 10 11 G (solution 2), respectively. Since 
such low E pulsars are generally very stable rotators with 
little timing noise or glitch activity, it was deemed likely 
that one of these is the true solution, and the other one 
is an alias with an incorrect cycle count. It is also im- 
portant that no solutions with smaller dipole B s , and no 
spinning-up solutions, were found. 
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TABLE 8 

XMM-Newton Phase- Resolved Spectra of PSR J0821-4300 



Model 


Group 1 


Group 2 


Group 3 


Group 4 Sum 


Parameter 


2001 Apr-Nov 


2009 Dec - 2010 May 


2010 Oct-Dec 


2011 Apr - 2012 Apr 2001 Apr - 2012 Apr 



kT w (keV) 
kT h (keV) 
F[A0i] a 
F[A<f> 2 ] a 
E (keV) 
a (eV) 
EW[k<t>i] 
EW[A<j> 2 ] 
xS(DoF) 



(eV) 
(eV) 



0.23(0.21,0.25) 
0.41(0.39,0.44) 
4.00(3.97,4.08) 
4.33(4.28,4.38) 
0.79(0.76,0.81) 
< 62 
40(24, 51) 
14(3,22) 
1.01(87) 



0.25(0.24,0.26) 
0.45(0.43,0.47) 
3.91(3.88,3.93) 
4.37(4.35,4.40) 
0.71(0.68,0.73) 
44(21,70) 
77(59,89) 
34(23, 42) 
1.19(229) 



0.24(0.23,0.25) 
0.44(0.42,0.46) 
3.94(3.91,3.96) 
4.37(4.34,4.39) 
0.72(0.69,0.74) 
68(45, 95) 
58(46, 76) 
42(32, 57) 
1.05(269) 



0.27(0.25,0.28) 
0.47(0.44,0.52) 
3.99(3.96,4.01) 
4.45(4.42,4.47) 
0.69(0.63,0.73) 
133(78, 195) 
86(67, 107) 
39(26,52) 
0.99(250) 



0.25(0.24,0.25) 
0.44(0.43,0.45) 
3.96(3.94,3.97) 
4.40(4.38,4.41) 
0.72(0.70,0.73) 
69(52,89) 
61(53,74) 
31(26,39) 
1.20(212) 



Note. — Results from simultaneous fits to XMM-Newton spectra of PSR J0821— 4300 extracted from two pulse-phase intervals, 
A(/>i, A<^>2, as defined in the Figure [5] Groups numbers arc defined in Tabic ]7] The blackbody temperatures and Gaussian line 
energy are linked between the two phases. The column density is fixed at A^h — 4.28 X 10 21 cm -2 , the phasc-avcraged value for the 
summed spectrum in Table [6] Quoted uncertainties arc la for three interesting parameters. 

a Absorbed flux quoted for the 0.5 — 5.0 keV band in units of 10 -12 erg s^ 1 cm~ 2 . 
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PSR J1 21 0-5226 




PSR J1210-5226 



0.5 - 2.5 keV 



A CHANDRA ACIS CC-MODE 
A XMM EPIC PN 



1000 2000 3000 4000 5000 

TIME (DAYS SINCE MJD 51549) 

Fig. 12. — Pulse-phase residuals from the linear term (dash-dot 
lin e) o f the phase ephemeris of PSR J1210— 5226 presented in Ta- 
ble [H Included are the new observations liste d in Table [9] and the 
archival timing observations from Table 1 of Halpcrn & Gotthclf 
l|201H ). The quadratic term (solid line) corresponds to the uniquely 
determined period derivative spanning the years 2000-2012. The 
error bars are generally smaller than the symbol size. 



2.2*10* - 




PULSE PHASE 

Fig. 13.— Pulse profiles of PSR J1210-5226 in the 0.5 - 2.5 keV 
band using data from a ll ti ming observations, folded according to 
the ephemeris in Table 1101 Phase zero corresponds to the listed 
TDB epoch of the ephemeris. Included are the new observations 
in Table [9] and the archiv al timing observations from Table 1 of 
IHalpern & GottheU 12011]) . 

these two CCOs both have velocities toward the high end 
of the distribution of pulsars, it would be a tantalizing 
physical result in itself. However, examination of another 
CCO, the NS in Cas A, doesn't neces sarily support such 
high v eloci ties for CCOs, in general; iThorstensen et all 
pOOll) and lFesen etHI ([20061 ) find v± « 350 km s" 1 with 
respect to the explosion center of Cas A. 



In 2011-2012 we started a new series of observations of 
PSR J1210-5226 with Chandra and XMM-Newton that 
was designed to initiate and maintain a unique, phase- 
connect timing solution of this pulsar for the first time 
and eliminate the prior timing ambiguity. Table [9] is a 
log of the new observations. The instrumental setups 
and analysis methods are the same as those described in 
Section 3 for PSR J0821-4300, except that the XMM- 
Newton source photons were extracted from a 20" radius 
aperture instead of 30", and the 0.5 — 2.5 keV band was 
selected to optimize the pulsed signal. 

The coherently measured period from the new ob- 
servations is precise enough to reject previous timing 
solution 2, while solution 1 with the smaller B s ex- 
trapolates precisely to the new period. Furthermore, 
the new pulse phases are aligned accurately with solu- 
tion 1, while solution 2 gives random phases. Thus a 
unique quadratic ephemeris fits all 12.9 years of Chan- 
dra and XMM-Newton timing. Table ITUl gives this global 
ephemeris, derived from a least-square fit to the TOAs 
as described in Section 3, with the phase residuals shown 
in Figure 1121 The previous uncertainty on the / of so- 
lution 1 is reduced by a factor of 2, with the result 
/ = -1.2363(91) x 10" 16 Hz s -1 ; the corresponding 
dipole magnetic field is B s = 9.8 x 10 10 G. Figure QJ 
shows the pulse profile using data from all the observa- 
tions folded on the presented ephemeris. 

Strictly speaking, the derived period derivative is an 
upper limit to the intrinsic one, as any proper mo- 
tion has not been measured and taken into account. 
PSR J1210-5226 and PSR J0821-4300 are at similar 
distances, and the location of PSR J1210— 5226 with re- 
spect to the geometric center of SNR PKS 1209-51/52 
allows (but does not require, since the kinematics of the 
SNR are unknown) a proper motion of ~ 70 mas yr -1 , or 
v±_ ~ 730 km s _1 (Dc Luca et all 120111 ). similar to that 
of PSR J0821-4300. If so, the Shklovskii effect given 
by equation (1) would contribute i\ ~ 1.1 x 10 -17 , or 
half of the total P, and the spin-down dipole field would 
be reduced somewhat, to B s w 7 x 10 10 G. Of course, if 



6. DISCUSSION 

6.1. On the Age of Puppis A 

Beck er et al.l (|2012t ) discussed the impact of the revised 
proper motion of PSR J0821— 4300 on the inferred age of 
Puppis A. The age of the SNR had been derived previ- 
ously as 3700 ± 300 yr from the proper motions of optical 
filaments that point back to a common center, presumed 
to be the site of the explosion, and assuming no decel- 
eration. The motion of the NS also extrapolates to the 
same center, but the distance traveled, 371" ± 31", cor- 
responds to an age of 5200 ± 1000 yr (or 6 100 ± 1000 yr 
for ou r proper motion of 61 ± 9 mas yr -1 ). iBecker et al.l 
(|2012t ) refer to these marginally contradictory measure- 
ments as independent, and choose to average them, giv- 
ing the value 4.5 kyr listed in Table [T] However, they are 
not truly independent, because they assume the same 
starting location. Furthermore, the discrepancy is worse 
if the filaments have decelerated. Just as the new Chan- 
dra observations of the NS have improved the accuracy 
of its proper motion, we suggest that a contemporary 
observation of the optical filaments of Puppis A may im - 
prove the precision of the orig inal I Winkler eTaD (1988) 
study on which the optical proper-motion age is based. 
Such an investigation could lead to a more detailed un- 
derstanding of dynamics of the filaments, and a more 
accurate age for Puppis A. 

6.2. PSR J0821-4300 and PSR J1210-5226 as 
Anti-Magnetars 

The spin-down of power PSR J0821-4300, E = 1.9 x 
10 32 erg s" 1 , is consistent with being caused by mag- 
netic braking of an isolated neutron star with a weak 
dipole field of B s — 2.9 x 10 10 G. Its spin-down power 
is much smaller than its observed thermal X-ray lumi- 
nosity, L x 5.6 x 10 33 d\ 2 erg s _1 , which rules out 
rotation as a signifi c ant po wer source. As discussed in 
IHaloern fc GottheU (|2010aD for PSR J1852+0040, the 
very small P disfavors propeller spin-down and accre- 
tion as a source of the X-rays. Thus, residual cooling 
remains the most plausible source of the X-ray luminos- 
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TABLE 9 

Log of New X-ray Timing Observations of PSR J1210— 5226 



Mission 


Instr/Mode 


UbslD 


Date 


Exposure 


btart rLpocn 


Frcqucncy a 


72 
Z l 








CUT) 


(ks) 


(MJD) 


(s) 


Chandra 


ACIS-S3/CC 


14199 


2011 Nov 25 


31.0 


55890.233 


2.3577625(28) 


48.2 


Chandra 


ACIS-S3/CC 


14202 


2012 Apr 10 


33.0 


56027.637 


2.3577709(43) 


22.2 


XMM 


EPIC-pn/SW 


0679590101 


2012 Jun 22 


26.5 


56100.537 


2.3577687(25) 


68.3 


XMM 


EPIC-pn/SW 


0679590201 


2012 Jun 24 


22.3 


56102.752 


2.3577621(34) 


69.0 


XMM 


EPIC-pn/SW 


0679590301 


2012 Jun 28 


24.9 


56106.490 


2.3577636(23) 


109.6 


XMM 


EPIC-pn/SW 


0679590401 


2012 Jul 02 


24.5 


56110.918 


2.3577626(28) 


61.4 


XMM 


EPIC-pn/SW 


0679590501 


2012 Jul 18 


27.3 


56126.553 


2.3577640(27) 


51.9 


XMM 


EPIC-pn/SW 


0679590601 


2012 Aug 11 


27.3 


56150.408 


2.3577637(23) 


81.8 


Chandra 


ACIS-S3/CC 


14200 


2012 Dec 01 


31.1 


56262.095 


2.3577634(28) 


39.8 



3aryccntric frequency derived from a Z x test. The given uncertainty is for the la confidence interval. 



TABLE 10 
Ephemeris of PSR J1210- 



5226 



Parameter 



Value 



R.A. (J2000) a 
Decl. (J2000) a 

Ephemeris Epoch (MJD TDB) b 
Ephemeris Span (MjD) 
Frequency, / 
Frequency derivative, / 
Period, P 

Period derivative, P 
Surface dipole dipole field, B 3 
Spin-down luminosity, E 
Characteristic age, t c 



12 n 10 m 00 s .91 
-52°26'28".4 
53562.0000006 
51,549-56,262 
2.357763502865(65) Hz 
(-1.2363 ± 0.0091) X 10" 16 Hz s 
0.424130748816(12) s 
(2.224 ±0.016) x 10" 17 
9.8 x 10 10 G 
1.2 X 10 31 erg s" 1 
302 Myr 



a Measured from Chandra ACIS-S3 ObsID 3913. Typical un- 
certainty is 0. 6. 



Epoch of fitted minima of summed pulse profile; phase zero 
in Figure [13] 

ity of PSR J0821-4300 and, by extension, of all CCOs. 
Given the meager spin-down power of PS R J0821— 4300, 
we ca n als o discount the suggest ion of iRevnoso et al.l 
( 2003) and lCastelletti et al.l (|2006l) that structure in the 
Puppis A SNR and associated H I is caused by jets emit- 
ted by the pulsar, similar to the SS433/W50 system. Par- 
enthetically, if the H I morphology is not caused by the 
pulsar, then it does not provide supporting evidence of 
its distance. 

An unresolved question about PSR J0821— 4300 is the 
origin of its phase-dependent, possibly variable emission 
or absorption feature. In the absence of accretion, it 
would be difficult to understand how an emission line 
is generated, or why it would vary in a few years. The 
indication of variability is not strong, and there is no 
other evidence of accretion. Therefore, we will consider 
here an absorption-line interpretation, for which there 
is a good precedent in PSR J1210— 5226, the only iso- 
lated pulsar t o show a series of strong absorption lines in 
its spectrum dSanwal et al.ll2002t iMereghetti et~aT]|2002l: 
iBignami et al.ll2003D . 

The spectral features in PSR J1210— 5226 are now 
widely considered to comprise the electron cyclotron fun- 
damental at Eq = 0.7 keV and its harmonics. The sur- 
face magnetic field strength where the lines are formed 
is inferred to be B 8 x 10 10 G according to equa- 
tion (2) and assuming z 0.3. The relative strength 
of the harmonics is explained by treating them as reso- 



nances in the photospheri c free- free opacity in the pres - 
ence of the magnetic field (jSuleimanov et al.1l2010L[20T^) . 
PSR J1210— 5226 is exceptional in having the largest 

^Kiown spin-down dipole magnetic field among CCOs, 

-4iew confirmed as B s < 9.8 x 10 10 G. This is only slightly 
larger than Bk8x 10 10 G inferred from the spectral fea- 
tures, and is the first case in which such a comparison has 
been made between independent methods of measuring 
surface B-field on an isolated pulsar. Effects that could 

-eliminate the already minor discrepancy include the un- 
measured proper motion, which would decrease the in- 
ferred spin-down dipole field, and the results of numerical 
models of a force-free magnetosphere, which imply a dif- 
ferent spin-down law from the standard vacuum dipole 

-expression 



Ze'IPP 



1/2 



in 2 R 6 sin 2 a J 
such that, more accurately, 



3.2 x 10 1 



PP 
sin 2 a 



1/2 



G 

(3) 



1/2 



47T 2 i? 6 [1 



sin 2 a| 



= 2.6x10 



id 



PP 



1 



sin 2 Q! 



(4) 

(|Snitkovskvl[27)0rl . This means that a measured spin- 
down rate is obtained with B s smaller by a factor 0.58 — 
0.82 than the conventional approximation B s = 3.2 x 

10l9(PP)l/2 Q 

In the case of PSR J0821— 4300, by the same logic, its 
weaker inferred dipole magnetic field of B s < 2.9 x 10 10 G 
wouldn't naturally account for a spectral feature at 
0.7 - 0.8 keV, instead predicting E < 0.26 keV. The 
absorption model with the cyclotron fundamental at 
Eq = 0.46 keV comes closer to the prediction. How- 
ever, additional variables can affect the comparison of 
a dipole field inferred from spin-down with that from 
a cyclotron line. First is the factor of 2 variation of 
dipole field strength over the surface, with B p = 2B S . 
Second, the actual field is not likely to be a centered 
dipole, and may be larger than B s or B p in places where 
the line is formed. The asymmetric distribution of sur- 
face temperature on PSR J0821— 4300 already appears 
to require a complex magnetic field geometry, such as 
an off-center dipole or higher multipoles. Third, the in- 
ferred spin-down field depends on the uncertain NS mass 
and radius as in equations (3) and (4). B s scales as 
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cx M 1 ' 2 R 



2 {l + z) (IRavenhall fc Pethick|[l99l . 

{1-2GM/RC 2 )- 1 / 2 cx B/E from equation 



uansporu in une crusi anu envelo pe 01 j 
modelled, most rec e ntly bylGeppert et alJ 
IPerez-Azorm et aT] (|2006al ) . and iPons et ; 



I 1 / 2 R- 
where 1+z ■ 

(2). For an astrophysically likely M — 1.4 M Q , theoreti- 
cal NS equations of state allow 8 < R < 15 km, therefore 
1.18 < 1 + z < 1.44. The spectroscopic B is there- 
fore uncertain by ±10% when we adopt z = 0.3, while 
the dipole B s is uncertain by the much larger factor of 
~ 2. (The gravitational redshift cancels out in their ra- 
tio.) Allowing for these uncertainties, it is reasonable to 
adopt the hypothesis that feature(s) in the spectrum of 
PSR J0821— 4300 are due to the cyclotron process. 

The continuum X-ray spectrum and pulse profiles of 
PSR J0821— 4300 are indicative of antipodal hot spots 
of dif f erent areas and temp eratures (jGotthelf fc Halpera 
20091 IGotthelf et all [20101. which are difficult to ac- 
count for if the magnetic field is weak. A related 
problem is t he high pulsed fraction of 64% from PSR 
J1852+0040 (IHalpern fe GottheU l2010ah . a timing twin 
of PSR J0821-4300. Polar cap heating by any mag- 
netospheric accelerator must be negligible as a source 
of surface heating in CCOs, being only a small frac- 
tion of the already insignificant spin-down power. Ther- 
mal X-rays from residual cooling can be nonuniform if 
there is anisotropic heat conduction in the star. The 
effect of different magnetic field configurations on heat 
transport in the crust an d envelope of NSs has bee n 

" J200U2006), 
all (|2009[) . A 

toroidal field is expected to be the initial configuration 
generated by differential rotation in the proto-neutron 
star dynamo ([Thomps on fe Duncan! 119931 ) . One of the 
effects of crustal toroidal field is to insulate the magnetic 
equator from heat conduction, resulting in warm spots 
at the poles. The warm regions can even be of different 
sizes due to the antisymmetry of t he poloidal compo- 
nent of the field ([Geppert et al.|[2006l ). which is evocative 
of the antipodal thermal structure of PSR J0821-4300. 
To have a significant effect on the heat transport, the 
crustal toroidal field strength required in these models 
is ~ 10 15 G, many orders of magnitude greater than 
the poloidal field if the latter is measured by the spin- 
down. Purely toroidal or poloidal fields are thought 
to be unstable in an initi ally fluid NS ([Tavlerl U973I: 
iFlowers fe Rudermanl [l977l) . although the toroidal field 
may be stabilized by a poloidal field tha t is several orders 
of magnitude weaker (jBraithwaite 2009). We suggest the 
la tter as a viable c o nfigu ration for a CCO. 

iShabaltas fc Lail ([2012[ ) tried to model the pulse pro- 
file of PSR J1852+0040 with anisotropic conduction, and 
concluded that they needed a toroidal crustal field of 
B<p > 2 x 10 14 G to produce its high pulsed fraction. 
Even then, they observe, the shape of the modelled light 
curve doesn't match the observed one. Since the shape 
of the light curve is not reproduced, the surface temper- 
ature distribution of PSR J1852+0040 and its physical 
o rigin are still not known. 

IPage fc Sarmientol (1996D IPerez-Azorm et all (|2006b[) . 
IZane fc Turollal (|2006l ). and lZand (|2007l ) investigated NS 
surface emission patterns using a combination of star- 
centered dipole and quadrupole magnetic field compo- 
nents to model asymmetric pulse profiles. The behavior 
of PSR J0821— 4300 may ultimately be explained by simi- 
lar models. It remains to be shown that if CCOs can have 



field configurations that are strong enough to affect heat 
transport to the extent required, while not exceeding the 
spin-down limits on the external dipole field. 

6.3. Origin and Evolution of Anti-magnetars 

PSR J0821-4300 is nearly a twin of PSR J1852+0040 
in its spin properties, and there are no other young 
neutron stars with measured magnetic fields this weak. 
They fall in a region of the P — P diagram (Fig- 
ure [3J that is devoid of ordinary (non-recycled) ra- 
dio pulsars (Manchester ct al. 2005|). They overlap 
with the supposed mild ly recycled pulsars in this area 
(jBelczvnski et all [2010). The characteristic age of 
254 Myr for PSR J0821-4300 is not meaningful because 
the pulsar was born spinning at its current period, as 
is the case for the other CCO pulsars. If their magnetic 
fields remain constant, they will not move in P — P space 
for > 10 s years (but see below). That CCOs are found 
in SNRs in comparable numbers to other classes of NSs 
implies that they must represent a significant fraction of 
NS births. If PSR J0821-4300 and PSR J1852+0040 are 
typical CCOs, the area around them in the P—P diagram 
should be densely populated with "orphaned CCOs" that 
remain after their SNRs dissipate in ~ 10 5 years. Why 
there are few ordinary radio pulsars and no older X-ray 
pulsars n ear the i r loca tion is then a mystery, as empha- 
sized by [Kaspi (20l"o|). which may indicate that radio 
luminosity is a function of spin-down power. 

There are not yet enough CCOs to know whether 
they are intrinsically radio-quiet relative to ordinary and 
recycled pulsars, rather than unfavorably beamed. It 
is also possible that some ordinary radio pulsars with 
B s < 10 11 G are actually much younger than their 
characteristic ages, and may be orphaned CCOs that 
could be recognized in X-rays. Whether or not they 
are radio pulsars, nearby orphaned CCOs should be de- 
tectable as thermal X-ray sources for 10 5 — 10 6 years, sim- 
ilar to the se ven ROSAT discovered isolated NSs (INSs: 
Habcrl 200 3) which ; however, have strong magnetic fields 
(jKaplan et all 120091) . It is likely that the INSs are kept 
hot for longer than CC Os by continuing magnetic field 
decay (|Pons et al.ll2007D . which would explain their ob- 
served abundance relative to the elusive orphaned CCOs. 
It may be difficult to detect and/or recognize orphaned 
CCOs if they cool faster than ordinary neutron stars. 
One effect that can accelerate cooling is an accreted 
light-element envelope, wh ich has higher heat c onduc- 
tion than an iron surface (jKaminker et all 120061 ). The 
newly discovered 59 ms pulsar 1RXS J141256.0+792204 
("Calvera"), originally detected in the ROSAT All-Sky 
Survey, may b e the first recogni z ed example of an or- 
phaned CCO (jZane et all [20lTl iRutledge et "all [200l 
iHalpernl 120111 ). pending a measurement of its period 
derivative. The isolate d NS 2XMM J10 4608.7-594306 
in the Carina Nebula (|Pires et al.l I2012T) has been sug- 
gested as another orphaned CCO. 

It is possible that the weak magnetic field of CCOs is 
causally related to slow rotation at birth thro ugh the tur- 
bulent dynamo ([Thompson fc Duncan! 1 19931 ) that gener- 
ates the magnetic field. (See Spruit 2008 for a review of 
possible mechanisms for the origin of magnetic fields in 
neutron stars.) If the dipole magnetic field B s is simply 
related to the initial spin period, and assuming that the 
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present period P is in fact the birth period because the 
spin-down time is so much longer than the true age, we 
may expect an anti-correlation between B 8 and P. With 
only three data points to compare, and with two of them 
having nearly identical values, there is not much evidence 
to examine for a trend. In fact, the CCO with the longest 
period, PSR J1210— 5226, also has the strongest dipolc 
field of the three, which would not by itself support such 
a simpl e correlation. A population analy sis of radio pul- 
sars by iFaucher-Giguere fe Kaspil (|2006l) concludes that 
there is a wide distribution of birth periods, with a mean 
of ~ 300 ms and a dispersion of a ~ 150 ms. If this is 
true, the birth periods of CCOs are not in fact long, and 
their weak dipole fields may be the effect of some as-yet 
unknown parameter. 

In an alternative theory for CCOs, a normal (~ 10 12 G) 
magnetic field is buried in the core or crust of a NS by 
prompt fall-back of supernova debris, and takes thou- 
sands of years to diffuse back to the surface, during 
which time the NS appears as an anti-magnetar. This 
is assuming that the accreted matter is itself not mag- 
netized. The timescale for diffusion is highly depen- 
dent on the amount of matter accreted. According to 
IMuslimov fc Page) (|1995D . for accretion of ~ 1O" 5 M , 
the regrowth of the surface field is largely complete after 
~ 10 3 yr, but if > 0.01 M© is accrete d, then th e diffu - 
sion time could be millions of years. IChevalierl (p89j) 
calculated that the neutron star in SN 1987A could have 
accreted ~ 0.1 Mq of fallback material in the hours after 
the SN explosion, aided by a reverse shock from the he- 
lium layer of the progenitor. If so, it may never emerge 
as a radio pulsar. 

Interesting support for the theory of field burial and 
regrowth is the absence of evidence for magnetic field 
strengths < 10 11 G in accreting high-mass X-ray binary 
pulsars, as inferred from their puls e periods and period 
derivatives ()Popov fc Turollall2013f) . If intrinsic fields in 
the range 10 10 — 10 11 G were common, equilibrium spin 
periods of 0.1 — 1 s should be frequent in HMXBs, but 
they are not. If the NSs in these systems are born in 
the same way as isolated pulsars, this could imply that 
birth fields are never so small; instead, field regrowth has 
occurred in all cases after it was initially buried. One 
caveat here is that ~ 40% of HMXBs do not have mea- 
sured spin periods, although there is no strong selection 
e ffect against detect ing P < 1 s. 

iBernal et al.l (|2010D has revisited the process of hyper- 
crit i cal ac cretion onto a magnetized neutron star, while 
Hoi ([20H made new calculations of the subsequent dif- 
fusion to constrain the magnetic fields of CCOs at birth 
and the accreted mass, finding 10~ 4 — 10~ 5 M Q for the 
latter. These models are difficult to test using CCOs, be- 
cause it would involve measuring the braking i ndex or the 
change of the dipole magnetic field directly. iHol ([201 1[ ) 
implied that the isolation of CCOs in the P — P diagram 
may be less problematic in this model because, as their 
dipole magnetic fields increase, they will evolve to join 
the bulk of the pulsar population. However, for the first 
~ 10 5 yr, field growth can only move a CCO vertically 
upward in the P — P diagram, as the braking index has a 
large, negative value. Orphaned CCOs should still have 
periods of ~ 0.1 s and lie in a sparse region. For this 
reason, an important test for an example of field growth 



is to determine the P of Calvera, which will reveal if its 
dipole field is greater than that of PSR J0821-4300 and 
PSR J1852+0040. 

Previous calculations of field burial and re-emergence 
for CCOs were one dimensional. The first two- 
dimensional calculatwns_Jor this p urpose were recently 
reported bv lVigano fc Pons! \20m . They find that the 
accretion must be essentially isotropic to bury the field 
sufficiently to cause the required orders-of-magnitude re- 
duction of the external dipole. If the accretion were in- 
stead confined to the equator or the magnetic poles, the 
reduction of the dipole component would not be signifi- 
cant enough to produce a CCO. Finally, it remains to be 
seen if the resulting surface thermal distribution during 
the CCO phase can be m ade compatible with o bserved 
spectra and pulse profiles. iVigano fc Pons! ([2012h briefly 
presented the temperature distribution from one of their 
models that produces hot polar caps at a time when the 
external dipole field is 10 10 G. But even in this case, the 
temperature does not vary by as much as a factor of 2 
over the surface, a range that would be required to match 
the properties of PSR J0821-4300 and other CCOs. 

7. CONCLUSIONS AND FUTURE WORK 

Measurement of the spin-down rate of the 112 ms 
PSR J0821-4300 in Puppis A was achieved using X- 
ray observations coordinated between XMM-Newton and 
Chandra, the resulting phase-connected ephemeris span- 
ning 2.3 yr. We also measured the proper motion of the 
pulsar in Chandra HRC images over 10.6 years. The 
proper motion makes a non-negligible contribution to 
the period derivative via the Shklovskii effect, and the 
uncertainty in proper motion and distance limits the 
accuracy of the intrinsic period derivative to rs 16%. 
The straightforward interpretation of these results is 
dipole spin-down due to a weak surface magnetic field 
of B s = 2.9 x 10 10 G, the smallest measured for any 
young neutron star, and nearly identical to that of the 
CCO pulsar PSR J1852+0040 in Kes 79. 

A phase-dependent spectral feature in 
PSR J0821— 4300 can be modelled either as an emission 
line of energy w 0.75 keV, or as a cyclotron absorption 
line and its harmonic, with Eq s» 0.46 keV. For reasons 
that are not clear, it is stronger during the pulse-phase 
interval in which the continuum spectrum is softer. 
There is only marginal evidence for long-term variability 
of this feature. The local magnetic field strength in the 
area where the spectral feature is produced may have to 
be larger than the dipole spin-down value. 

The existing spin-down measurements for three CCOs, 
including a definitive new result for PSR J1210— 5226, 
are compelling evidence that a weak dipole field com- 
ponent is the physical origin of the CCO class in gen- 
eral. It is reasonable to assume that the CCOs which 
have not yet been seen to pulse have magnetic fields that 
are similar to or weaker than those of PSR J0821— 4300 
and PSR J1852+0040. Otherwise, if their fields were 
stronger, their spectra should show cyclotron lines simi- 
lar to PSR J1210-5226 and (possibly) PSR J0821-4300. 
Deep X-ray timing searches of the remaining members of 
this class, to smaller limits on pulsed fraction, could still 
discover new pulsars and supply valuable data on their 
birth properties and evolution. 

A remaining theoretical puzzle about CCOs is the ori- 
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gin of their surface temperature anisotropics, in partic- 
ular, the one or two warm/hot regions that are smaller 
than the full neutron star surface. The measured spin- 
down power is too small to contribute to this emis- 
sion. Continuing accretion is unlikely because of the 
small spin-down rates. It appears that any explana- 
tion will have to involve stronger magnetic field com- 
ponents in the crust, with toroidal or quadrupolar ge- 
ometries, that do not contribute to the dipole spin-down 
torque. A physical model of thermal and magnetic field 
structure that self-consistently reproduces the spin-down 
rates, X-ray spectra, and pulse profiles, is still needed. If 
it can be studied with better data, the details of the 
phase-dependent and possibly variable spectral feature 
from PSR J0821— 4300 may contribute to a more defi- 
nite model of its surface magnetic field geometry. 

The CCO pulsars PSR J1852+0040 and 
PSR J082 1-4300 fall in a region of B - P space 
that overlaps with what are assumed to be moder- 
ately recycled pulsars, but is otherwise empty. An 
understanding of the evolutionary status of CCOs is 
critically dependent upon a search for their descendants, 
the orphaned CCOs without SNRs. Some single radio 
pulsars with spin parameters similar to CCOs may be 
orphaned CCOs rather than recycled pulsars, and they 
may be much younger than their characteristic ages. 
X-ray observations of pulsars in this region may find 
evidence of their relative youth via surface thermal 



emission. The spin-down rates of orphaned CCOs would 
reveal whether CCOs have intrinsically weak magnetic 
fields, or if field re-emerges on a timescale of ~ 10 4 yr 
from a normal magnetic field that was buried by prompt 
fallback of supernova debris, as is invoked in some 
theoretical studies. Such rapid evolution would lead to 
orphaned CCOs that lie directly above CCOs on the 
P — P diagram but are still detectable as thermal X-ray 
sources. The radio-quiet pulsar Calvera is a possible 
such candidate. Whether or not CCOs have buried 
fields, the paucity of radio pulsars with similar spin 
parameters is real and requires an explanation. 
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